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The Pierre Auger Observatory infers the chemical composition of ultra-high-energy cosmic 
rays through two independent detection techniques. The Fluorescence Detector (FD) measures 
the longitudinal profile of high energy air showers and can determine the depth of the shower 
maximum -Y max , which is sensitive to the chemical composition of the primary cosmic rays. 
Additionally, measurements by the Surface Detector (SD) provide independent experimental 
observables based on the muonic shower component to analyze the chemical composition. 

We present the results for the A' max distributions and the mass composition results measured 
by the FD and the SD for the energies E > 10 18 eV. The data will be compared with the 
expectations for proton and iron primaries according to different hadronic interaction models. 


1 Introduction 

The chemical composition of ultra-high-energy cosmic rays (UHECRs) with energies E > 10 18 eV 
can only be studied indirectly by measuring the secondary particles produced in extensive air 
showers. The determination of the composition is important to understand the origin of UHE¬ 
CRs. In this work, we show some of the latest measurements of composition related observables 
of the Pierre Auger Observatory in comparison with different interaction models. 

1.1 Pierre Auger Observatory 

The Pierre Auger Observatory is the world largest detector for UHECRs with an aperture of 
3000 km 2 covered by the Surface Detector array 1 (SD) of 1660 water-Cherenkov stations and 
overlooked by the Fluorescence Detector 2 (FD) with 27 telescopes at the borders of the array. 
This hybrid design with 2 independent detector components allows one to study simultaneously 
the longitudinal shower profile in the atmosphere and the lateral distribution of shower particles 
on the ground with high accuracy. Due to the low fluorescence light flux of extensive air showers 
the FD can only take data during clear, moonless nights, which corresponds to a duty cycle of 
13%, whereas the SD can measure with a duty cycle of nearly 100%. 

2 Studies of the Chemical Composition 

The Pierre Auger Observatory measures the primary particle composition with both detector 
components, the FD and the SD. In this work we present two different analyses. 
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Figure 1: Development of the mean depth of the shower maximum (A' max ) as a function of the energy, with a 

broken line fit of the elongation rate 5 . 




Figure 2: Energy evolution 5 of the average of the shower maximum (A' max ) (left) and event-by-event fluctuations 
cr(A max ) (right) compared to the high-energy interaction models Sibyll2.1 6 , EPOS-LHC 7 and QGSJetII-04 8 . 


2.1 Fluorescence Detector Analysis 

One way to determine the mass of the primary cosmic ray is to study the longitudinal shower 
profile of the electromagnetic component in the atmosphere with the FD and to measure the 
depth of the shower maximum A max , which is the position of the maximum of energy deposition 
per atmospheric slant depth of an extensive air shower 3,4 . Lighter primaries penetrate the 
atmosphere deeper than heavier primaries. Also, due to the larger number of nucleons and the 
larger cross section, the event-by-event fluctuations of X max should be smaller for heavier nuclei. 
So the first two moments of the X max distribution, which are the mean (A max ) and standard 
deviation cr(A max ) are used to study the mass composition in the energy range above 10 178 eV. 
In figure 1 the average (A max ) is compared with energy 5 . It shows that the elongation rate D, 
which is defined as 

jj _ d(X m3iX ) 

d log E ’ 

cannot be described by a single line fit ( x 2 /ndf = 128.1/16), but requires a broken line fit 
{x 2 /ndf = 10.3/14). This suggests a change in the composition of UHECRs with increasing 
energy. In figure 2 the energy evolution 5 of the average of the shower maximum (X max ) and 

















Figure 3: Reconstructed muon production depth (MPD) 10 of an extensive air shower with a zenith angle 9 = 
(59.06 ± 0.08)° and an energy E = (92 ± 3) EeV with a Gaisser-Hillas function fit 11 . 


event-by-event fluctuations <x(V max ) are shown in comparison with the recent interaction models 
for extensive air showers Sibyll2.1 6 , EPOS-LHC 7 and QGSJetII-04 8 . This comparison suggests 
that with increasing energy there is a transition from a lighter to a heavier mass composition 9 . 
However, this interpretation is strongly dependent on the interaction models. 


2.2 Surface Detector Analysis 

It is possible to study the primary composition of UHECRs by reconstruction of the muon 
production depth (MDP) along the shower axis 10 . Due to the fact that muons come from pion 
and kaon decays, the MPD allows to study the hadronic component of the shower. This can 
be inferred from the timing information of muons in the SD stations, and a geometrical back- 
projection to the shower axis. In this way, SD stations can be used to reconstruct the longitudinal 
shower profile. However, signals measured in the SD are a mix of the electromagnetic and the 
muonic component of the air shower. To reduce the electromagnetic contamination, only air 
showers in the zenith angle band 55° < 9 < 65° are included and to reduce distortions in the 
reconstructed depths, only SD stations which are in the range of 1700 m to 4000 m from the 
shower core are taken. By means of the longitudinal distribution of the MPD, the V^ ax , which 
is the maximum of produced muons per atmospheric slant depth, can be reconstructed by a 
fit with a Gaisser-Hillas function 11 , as shown in figure 3. Similar to the {V max ) study in the 
previous section, the mean (X^ ax ) based on Monte-Carlo studies is smaller for heavier and larger 
for lighter primary cosmic rays. In figure 4 the data of this study are shown for showers with 
an energy E > 20 EeV. The interpretation of these data in terms of chemical composition is 
very challenging due to the considerable differences in the muon production number 12,13 and 
elongation rate based on the interaction models as shown in figure 4. This study can be used to 
constrain the interaction models and improve the understanding of hadronic interactions. The 
slope of the elongation rate suggests a change in the composition of UHECRs. However, this 
study currently suffers from limited statistics and the resolution on the MPD. Therefore, the 
measurements currently cannot rule out a constant composition. 
















Figure 4: Energy evolution of (X r ^ lax ) compared with interaction models EPOS-LHC 7 and QGSJetII-04 8 . The 

gray rectangles represent the systematic uncertainties 10 . 


3 Future prospects 

The sensitivity of measurements to the chemical composition of UHECRs will increase consider¬ 
ably in the next few years through increased statistics and by improved analyses. Additionally, 
to improve the sensitivity to the composition, the Pierre Auger Collaboration plans to enhance 
the electronics of the SD stations to improve the MPD resolution, and to add a dedicated muon 
detection system to the ground array to perform an event-by-event measurement of the muonic 
longitudinal shower profile with with large statistics 14 . 
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